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We review the current status of analyses of multijet final states in e+e - annihi- 
lation. Results for jet observables from LEP 1, LEP 2 and from the reanalysis of 
the PETRA experiment JADE will be presented. A determination of the b-quark 
mass using jet observables will be discussed and tests of power correction models 
will be shown. Finally, determinations of the QCD colour factors from an analysis 
of event shape distributions at several energy points using power corrections will 
be discussed. 

1 Introduction 

Hadronic final states in e + e _ annihilation events are the subject of many 
experimental and theoretical studies. The structure of the hadronic final 
states in e + e _ annihilation is characterised by the presence of a small number 
of so-called jets, i.e. clearly separated and collimated sprays of particles. As a 
consequence, hadronic events may be classified by e.g. the number of jets after 
a jet finding algorithm has been defined. Alternative classification schemes 
are event shape observables, where the reconstructed momenta of the final 
state particles are combined in a way which characterises the structure of the 
event in a single number. 

Section |^ presents results on jet production for centre-of-mass (cms) en- 
ergies from 35 to 189 GeV, while section || describes a measurement of the 
mass of the b-quark at a scale of M^o . Section ^ contains a brief summary of 
experimental tests of power corrections. Section || shows a test of the gauge 
structure of QCD and, finally, section |^ gives a summary of the report. 

2 Jet Production from 35 to 189 GeV 

Jet production was studied comprehensively using data from the JADE and 
OPAL experiments at the PETRA and LEP e + e~ colliders at cms energies 
from 35 to 189 GeV |Q. The large range of cms energies allows detailed studies 
of scale dependent effects predicted by QCD. 

Figure [l] (left) shows the 3-jet fraction measured using the JADE algo- 
rithm H at cms energies 35 (triangles) , 91 (points) and 189 GeV (squares) 
as functions of the jet resolution parameter y cut . In the JADE algorithm the 
invariant masses my of all pairs of final state particles are calculated and 
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Figure 1. The figure on the left shows the 3-jet fraction using the JADE algorithm as 
a function of the jet resolution parameter y cu t at y/s = 35 (triangles), 91 (points) and 
189 GeV (squares). The lines display Monte Carlo predictions by PYTHIA (solid), HER- 
WIG (dashed), ARIADNE (dotted) and CO JETS (dash-dotted). The figure on the right 
displays the 3-jet fraction at y cu t = 0.08 corrected for experimental and hadronisation 
effects as a function of the cms energy y/s M . 



the pair with the smallest value is combined into a pseudo-particle by adding 
their 4-vectors. This process is repeated until no invariant masses below a 
cut value remain: mfj/s > y cu t for all i, j. The decrease of the 3-jet rate at 
large y cu t is clearly visible. The distributions are well described by Monte 
Carlo simulation programs tuned to OPAL data recorded at y/s — M z o |B0] 
including QCD coherence effects in the parton shower, PYTHIA || (solid), 
HERWIG 0] (dashed) and ARIADNE [§ (dotted). A simulation program 
without coherence effects, CO JETS § (dash-dotted), describes the data less 
well. 

Figure [l] (right) shows the dependence of the 3-jet fraction i?3(y cut = 0.08) 
using the JADE algorithm as a function of cms energy yfs, corrected for detec- 
tor and for hadronisation effects. In leading order, we have Rz{\fs) ~ a s (\fs) 
for y cu t = 0.08 jj); the data provide convincing evidence for the running of a s 
as required by QCD. 

3 Running b- Quark Mass 

QCD predictions are generally calculated for massless quarks. This is a good 
approximation for the light (O(10— 100) MeV) u-, d- and s-quarks. However, 
for the heavy c- and b-quarks the masses are comparable to energy scales 
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Figure 2. The figure on the left shows the ratio of 3-jet rates in b- and d-quark events cor- 
rected to the parton level compared with NLO QCD predictions and simulations. The fig 
on the right shows the final result for mt>(M z o) compared with other measurements 



where perturbative QCD calculations are expected to be valid and thus quark 
mass effects may be significant. Quark masses in the QCD Lagrangian are 
free parameters like the strong coupling a s and have to be renormalised to 
obtain finite predictions, see e.g. |l3f| . The renormalised quark masses are 
expected to "run" , i.e. to depend on the energy scale of the process, because 
they must obey a renormalisation group equation (RGE). The main effect of 
a heavy quark mass in QCD is the suppression of gluon radiation from the 
heavy quark which leads to the expectation of a reduced 3-jet rate. 

The ALEPH collaboration presented an analysis to test the prediction 
of a running b-quark mass based on jet observables determined for b- and 
light quark events |L4|]. After correcting for experimental effects, b-tagging 
biases and hadronisation effects, the ratio of the jet observable measurements 
in b- over d-quark events is calculated, e.g. for the 3-jet rate R3 one has 
(R3) — i?3,b/-R3,d- Figure H (left) presents J?{^ rt (i?3) at several values of 
y cu t compared with NLO QCD predictions for — 3 and 5 GeV and Monte 
Carlo simulations. The data clearly prefer the lower value of while the 
simulations are in slight disagreement at low y cu t- 

The final measurement of m\ 1 (Mzo) is performed using the observable 
with the smallest hadronisation corrections and systematic uncertainties; this 
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Figure 3. The figure on the left shows as solkLlines 0(cP s ) QCD fits with power corrections 
to fst moments of event shape observables [ |15| . The figure on the right shows a summary 
of all results for ao and a s (M z o) from various analyses [J15l4l9 . 



is the 1st moment of the differential 2-jet rate distribution. The result 
rab(Mzo) = 3.27 ± 0.52 GeV is presented in figure ^| (right) together with 
other measurements at the Z° scale and at low scales. The QCD prediction 
of a running b-quark mass starting from the value at low scale is in good 
agreement with the measurements at M-^o . Assuming the QCD description of 
heavy quark effects to be correct the ratio ag/a" dsc = 0.997 ± 0.009 is deter- 
mined and provides a precise test of the flavour independence of the strong 
coupling a s . 



4 Power Corrections 

Most measurements in QCD studies have to correct for the discrepancy be- 
tween perturbative QCD calculations and the quantities calculated from the 
observed hadrons. These corrections are commonly carried out using Monte 
Carlo models of the hadronisation process like JETSET/PYTHIA, HERW1G 
or ARIADNE. An alternative approach to the problem of hadronisation are 
analytical QCD based models of hadronisation, the power corrections. 

In the Ansatz of Dokshitzer, Marchesini and Webber (DMW) |po| , the 
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effects of gluons with transverse momentum k t ~ Aqcd, so-called gluers, are 
calculated. The model must assume that the strong coupling a s is finite in 
the region of the Landau Pole leading to a new free parameter ao in the 
model: ao — J^ 1 a s (k)dk. The variable /!/ is the infrared matching scale 
where non-perturbative and the perturbative evolution of a s are merged. 

For the differential distributions of the event shape observables Thrust, 
Heavy Jet Mass, C-parameter and Total and Wide Jet Broadening, the model 
predicts that hadronisation effects are described by a shift of the perturbative 
prediction: F(y) — Fp^(y—c y P) where y is the value of the observable pl| , [2^ ]. 
For the 1st and 2nd moment one obtains (y) = (y) PT + c y P and (y 2 ) = 
(y 2 )pT + 2(y) PT c a P + 0(1/Q 2 ). The quantity c y depends on the observable 
while P ~ At hi /Qopf/xj) is universal and the Milan factor At takes account 
of two-loop effects |21j| . The shift is inversely proportional to the hard scale 
Q usually identified with the cms energy. 

A study using 1st moments of event shape observables by DELPHI is 
shown in figure || (left) Jl5|] using DELPHI data from LEP 1 and 2 and data 
from various experiments at lower energies. The fits of 0{a 2 s ) QCD predic- 
tions with power corrections (solid lines) describe the data well. The dashed 
lines represent the perturbative part, such that it becomes apparent that 
hadronisation corrections are important even at large cms energies. 

A direct test of the power corrections using differential distributions of 
the event shape observables is presented in fTc| . Data measured at y/s = 
35 to 183 GeV are fitted simultaneously with only a s (M z o) and ao as free 
parameters. The fitted predictions describe the data well within the fitted 
regions. 

Results from power correction analyses for a.s{AIyp) and ao from many 
recent analyses are summarised in figure | (right) |§|9|^| The results for 
a s (M z o) are generally consistent with the world average value a s (A/ z o) = 
0.119 ± 0.003 jl2| while the results for ao are in agreement with each other 
at the 20% level, as expected theoretically ^lj. The results for B w from 
distributions are not as consistent with universality of ao as the other results. 



5 QCD Colour Factors 

A study of the QCD colour factors using fits of C(a^)+NLLA QCD predictions 
with power corrections to distributions of 1 — T, C, B>t and £?w was presented 
in [Q. The QCD colour factors rifT-p, Ca and Cf represent the relative 



a Results for ao based on the old erroneous value of the_Milan factor M = 1.795 have been 
scaled to correspond to the correct value of M = 1.49 [E3| 
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Figure 4. The figure on the left shows the results of fits to event shape distributions with 
a s (M z a), ctQ and one of the QCD colour factors nf, C\ or Cp as free parameters. The 
vertical dotted lines indicate the expectations from standard QCD. The figure on the right 
shog the combined result of simultaneous fits of ce s (M z o), olq, C\ and Cp to 1 — T and 



contributions of the QCD vertices of quark-pair production from a gluon, 
gluon-radiation of a gluon (triple gluon vertex) and gluon radiation of a quark, 
respectively. In the product rifTp, n? is the number of active quark flavours 
and Tp is the actual colour factor. The QCD colour factors are determined 
by the choice of the gauge symmetry group, SU(3) in the case of QCD, and 
are expected as nf = 5 for Tp = 1/2, Ca = 3 and Cf = 4/3. 

The analysis follows [|l8| but uses more recent data. The dependence of 
the complete QCD predictions on the colour factors is made explicit such 
that the colour factors rif, Ca or Cf can be varied in the fits, in addition to 
a s (Mzo) and ao. The main sensitivity comes from the running of a s , which in 
0(a s ) reads a s (Q) = a s (p)/(l - 2/3 a s (fj,) InQi/Q)), where (3 — IICa - 2n f . 
Using the power correction calculations as the hadronisation model reduces 
potential biases from hadronisation corrections, because the power corrections 
depend explicitly on the QCD colour factors. 

Figure [| (left) shows the results of fits with a s (M z o ) , ao and one of the 
colour factors rif, Ca or Cf as free parameters. The results for a s (Mza) are 
consistent with the world average, except the result from B-yv when fitting 
Cf • The results for ao are consistent with those presented in section |]. The 
results for the colour factors are consistent with the SU(3) QCD expectation 
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with five active flavours, indicated by the vertical dotted lines. 

Figure [| (right) displays the unweighted averages of results of 1 — T and 
C from simultaneous fits with a s (M z o), C'a and C-p as free parameters; the 
non-perturbative parameter ao was fixed at o.q = 0.543 ± 0.058. Also shown 
arc the expectations for various alternatives for the gauge symmetry group in 
QCD, in particular U(l) 3 is the representation for a theory with three different 
neutral gauge bosons in direct analogy to QED. The measurement agrees well 
with standard QCD with the SU(3) symmetry group. 

The analysis is complementary to the traditional approach of using angu- 
lar correlations in 4-jet final states |25|-p7| and has similar total uncertainties. 
Under the assumption that QCD based on the SU(3) gauge symmetry group 
is the correct theory of strong interactions, the analysis provides a successful 
consistency check of the power correction model. 

6 Summary 

We have shown experimental studies of jet production in e + c~ annihilation. 
Jet production as measured from PETRA to LEP 2 energies is well described 
by QCD models and by perturbative QCD calculations. A measurement of the 
b-quark mass at the Z° peak provided evidence for the running of the b-quark 
mass as predicted by QCD. Investigations of power corrections were discussed 
and it was found that the model successfully predicts the hadronisation effects 
for a number of event shape observables. The free non-perturbative parameter 
ao is observed to be universal within the theoretically expected uncertainty of 
about 20%. A measurement of the QCD colour factors using power correction 
calculations was presented. This analysis is complementary to traditional 
analyses of angular correlations in 4-jet final states at the Z° peak and of 
similar accuracy. 

The author would like to express his gratitude towards the organisers of 
this meeting for a stimulating conference in a pleasant atmosphere. 
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